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Abstract: 

We construct a representation of the Separated Variables (SoV) for the quantum SL(2, R) Heisen- 
berg closed spin chain following the Sklyanin's approach and obtain the integral representation 
for the eigenfunctions of the model. We calculate explicitly the Sklyanin measure defining the 
scalar product in the SoV representation and demonstrate that the language of Feynman di- 
agrams is extremely useful in establishing various properties of the model. The kernel of the 
unitary transformation to the SoV representation is described by the same "pyramid diagram" 
as appeared before in the SoV representation for the SL{2, C) spin magnet. We argue that this 
kernel is given by the product of the Baxter Q— operators projected onto a special reference state. 
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1. Introduction 



In this paper we study the spectral problem for quantum SL(2,M) spin magnet within the 
Quantum Inverse Scattering Method PQ 13 El- The model is a generalization of the famous 
spin— 1/2 Heisenberg XXX chain. It describes the nearest neighbour interaction between N spins 
with the corresponding spin operators being the generators of infinite-dimensional representation 
of the SL(2, R) group. 

The interest to studying noncompact spin magnets is twofold. From the one side, there 
exists a deep relation between integrable models and quantum (3 + 1)— dimensional Yang- Mills 
theories jl]. It turns out that the scale dependence of certain correlation functions on the light- 
cone, like those defining a baryon distribution amplitude and twist-three light-cone distributions, 
is governed by the evolution equations which upon redefinition of the variables coincide with the 
Schrodinger equation for the SL(2,M) spin chain Hj. The number of sites in the lattice 
model is equal to the number of fields entering the correlation functions. The spin operators 
are the generators of the so-called collinear SL(2, R) subgroup of the full conformal group of 
the classical Yang-Mills Lagrangian. Another example comes from the studies of high-energy 
(Regge) asymptotics of the scattering amplitudes in multi-colour QCD. As was shown in [HIE], 
the spectrum of multi-gluonic compound states responsible for a power rise with the energy of 
the scattering amplitudes is described by the SX(2,C) spin magnet. 

From the other side, an exact solution of the spectral problem for completely integrable 
quantum mechanical systems with infinitely dimensional quantum space, like periodic Toda chain 
and noncompact SL(2) spin magnet, represents a challenge for the theory of Integrable Models. 
The conventional methods like the Algebraic Bethe Ansatz PJ 0] (ABA) are not always applicable 
to such systems and one has to rely instead on the methods of the Baxter Q— operator JU] and 
the Separation of Variables (SoV) llj. Being combined together, the two methods allow us to 
determine the energy spectrum of the model in terms of the eigenvalues of the Q— operator and 
construct an integral representation for the corresponding eigenf unctions by going over to the 
representation of the separated variables. 

In spite of the fact that the methods have been formulated a long time ago, the explicit 
construction of the Q— operator and the unitary transformation to the SoV representation for 
each particular model remains an extremely nontrivial task (for some known examples see review 
[T2*] and references therein). It is equally nontrivial to solve the emerging functional relations 
(the Baxter equations) and reconstruct a fine structure of the spectrum. 

A powerful algebraic approach to constructing the SoV representation has been developed by 
Sklyanin [TT]. It allows one to establish the relations both for the eigenf unctions in the separated 
coordinates and for the scalar product in the SoV representation. Their solutions are defined up 
to multiplication by an arbitrary periodic function. For models with finite-dimensional quantum 
space, like conventional 577(2) Heisenberg magnet, the separated coordinates take a discrete, 
finite set of values and, as a consequence, the above relations can be uniquely solved. For 
models with infinite-dimensional quantum space, like SL(2,M.) Heisenberg magnet and periodic 
Toda chain, the separated coordinates take continuous values and the question arises how to 
fix the ambiguity or, equivalently, what are the additional conditions that one has to impose 
on the eigenf unctions and the integration measure in the SoV representation. To answer this 
question within the Sklyanin's approach, one has to provide an explicit construction of the unitary 
transformation to the SoV representation and identify the analytical properties of the resulting 
expressions for the eigenfunctions in the separated coordinates. 
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For the periodic Toda chain this program has been carried out by Kharchev and Lebedev in 
the series of papers [T3j . The approach proposed in Ref. [T3] is based on the relation between the 
wave functions of the periodic and open Toda chains originating from a specific form of the Lax 
operator. Notice that similar relation does not hold for the SL(2) magnet. In the present paper 
we construct the SoV representation for the SL(2, M.) spin chain following the method developed 
in [HI [THj in application to the SL(2, C) spin magnet. The method is general enough as it is 
applicable to both the SL(2,~R) spin chain and the Toda chain. 

Our main results include: calculation of the Sklyanin's measure defining the scalar product in 
the SoV representation, establishing the relation between the kernel of the unitary transformation 
to the SoV representation and the Baxter Q— operator constructed in [T"~], proof of the equivalence 
of the ABA and the SoV method for the SL(2, M.) spin chain. 

The paper is organized as follows. In Section 2 we define the model and describe its integrabil- 
ity properties. In Section 3 we apply the Sklyanin approach and present an explicit construction 
of the unitary transformation to the Separated Variables for the SL(2,M.) closed spin chain. We 
demonstrate that the kernel of the SoV transformation admits a simple interpretation in terms 
of Feynman diagrams. In Section 4 we follow the diagrammatical approach and derive the in- 
variant scalar product in the SoV representation. In Section 5 we establish the relation between 
the transition function to the SoV representation and the Baxter Q— operator. It allows us to 
obtain the expressions for the eigenf unctions of the model in the separated variables. Section 5 
contains concluding remarks. Some details of the calculations can be found in Appendix A. In 
Appendix B we argue that the Algebraic Bethe Ansatz and the SoV method lead to the same 
expressions for the eigenstates of the SL(2,~R) spin magnet. 

2. The quantum »SX(2, M) spin magnet 

The quantum 5L(2,R) spin magnet is a one-dimensional lattice model of N interacting spins 
S n = (5° , , S~) (with n — 1, N). The spin operators in different sites commute with each 
other and obey the standard si (2) commutation relations 

[S° n , 5f] = ±hSt , [S+, S-} =2hS° n . (2.1) 

The corresponding quadratic Casimir operator is defined as 

&n = Sn^ri + ^(^n^n + &n) = ^ S n( s n ~ 1) • (2.2) 

— * — # 

We shall impose periodic boundary conditions, Sn+i = Si, and put h = 1 for simplicity. In 
addition, we shall assume that the spin chain is homogenous, Si — ... — sjv = s. Generalization 
to the case of inhomogeneous spin chain will be discussed in Section 6. 

2.1. Hamiltonian of the model 

The Hamiltonian of the homogenous SL(2, M) spin magnet is defined as [2*1 IT7] 

N 

T-Cn = H n>n+1 , H n , n+ i = ijj{J n>n+ i) - ip(2s) , (2.3) 

n=l 
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where ip(z) = d\nT(z)/dz is the Euler ^-function. The operator J nyn+ i is related to the sum of 
two neighbouring spins 

Jn,n+l{Jn,n+l — 1) = {S n + S n+ i) , (2-4) 

Jn,n+i — Jn,i an d its eigenvalues satisfy the condition J ntn+ i > 1/2. The Hamiltonian ()2.3j) has 
been constructed in Ref. [2j as a generalization of the spin— 1/2 XXX Heisenberg spin chain to 
high-dimensional representations of the SU (2) group. By the construction, TCn possesses a set of 
mutually commuting integrals of motion that we shall denote as q = (q^, ...,gjv)- Their number 
is large enough for the model to be completely integrable. The definition of the operators q will 
be given below (see Eq. (|2.18|) ). 

In what follows, we shall use a particular representation for the spin operators 

— iz n p n -\- 2s z n , S n = — ip n , S n = iz n p n + s , (2-5) 

where p n = —id/dz n and [x n ,p m ] = iS nm . The spin s is assumed to be real and s > 1/2. 
In this representation, the SL(2, M.) spin magnet ()2.3|) can be interpreted as a one-dimensional 
quantum mechanical model of N interacting particles with the coordinates z n and the conjugated 
momenta p n (n = 1, N). At N = 3 this model has appeared in high-energy QCD as describing 
the spectrum of the anomalous dimensions of the baryon distribution amplitudes jSj. 

The spin operators ()2.5|) act on the Hilbert space, V n , of functions *f?(z n ) G V n holomorphic 
in the upper half-plane , Imz„ > 0, and normalizable with respect to the scalar product [TH] 

ll^ll 2 = / Vz\^(z)\\ (2.6) 

Jlmz>0 

with z = x + iy. Here integration is performed over the upper half-plane and the integration 
measure isdefined as 2 

Vz = ^Zld 2 z(2lmz) 2s - 2 = ^Zldxdy(2y) 2s - 2 . (2.7) 

7T 7T 

The spin s in Eqs. f!2.5|) and ()2.7|) takes arbitrary real values s > 1/2. For s integer or half integer, 
the Hilbert space V n coincides with the linear space of the unitary irreducible representation of 
the SL(2, R) group of the discrete series [TH] 

< 2 > - jdw- * (^) • » = ( W ) e Si < 2 - R > ■ 

In this paper we shall study the Schrodinger equation for the Hamiltonian (|2.3|) 

H N ^ q (zi, z N ) = E q z N ) . (2.8) 

The Hamiltonian Ti^ acts on the quantum space of the model Vn, given by the direct product 
of the Hilbert spaces in each site, Vn = Yln=i I n distinction with conventional (com- 

pact) Heisenberg spin magnet, Vat is infinite-dimensional for arbitrary finite N. The eigenstates 

1 One can choose instead $(z„) to be holomorphic in the lower half-plane. As we will show below, the two 
cases, Im z > and Im z < 0, correspond to the different values of the total momentum of the system, p > and 
p < 0, respectively. 

2 Performing the conformal mapping w = i(z — i)/(z + i), one can bring this expression to a canonical form 
involving the integration over an interior of the unit disk in the w— plane |18| . 
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"$?q(zx, Zn) G Vat are holomorphic functions of the z— coordinates in the upper half-plane, 
lmz n > for n — 1, N, normalizable with respect to the scalar product 

II^atII 2 = I V N z\V(z 1 ,...,z N )\ 2 . (2.9) 

Here j V N z = Yl n =i fi mz >o ^ Zn an< ^ ^ ne measure Dz n is given by ()2.7|) . In Eq. (|2.8|) we indicated 
explicitly the dependence on the integrals of motion q = (q 2 , gjv)- The Hamiltonian (|2.3|) 
commutes with the total spin of the magnet S = Si + ... + on an d one of its component, 
iS- = pi + ... +Pn, defines the total momentum of the system. This allows us to assign a definite 
value of the momentum to the solutions to ()2.8|) . iS^ q)P = p^ Q)P , leading to 

/oo 
dx e ipxo m q ( Zl - x , ...,z N - x ) , (2.10) 
-oo 

where integration goes along the real x— axis. Notice that in virtue of the SL(2) invariance of 
the Hamiltonian, [Hn, S] = 0, the energy E q does not depend on p. 

It is straightforward to verify that for arbitrary real s the spin operators ()2.5j) are anti- 
hermitian with respect to the scalar product ()2.(i|l and (|2.9|) 

This property ensures that the Hamiltonian ()2.3j) and the total momentum operator iS- are 
hermitian on the Hilbert space of the model. As a consequence, the energy E q and the total 
momentum p take real values. As we will show in Section 3, the property (|2.11|) plays an 
important role in our construction of the SoV representation. 



2.2. Integrals of motion 

To construct the integrals of motion of the model, one follows the R— matrix approach QjJ |SJ |2J- 
The Lax operator for the SL(2,M.) magnet is defined as 

L n {u) =u + z(a.S n )=( U + f ° u ) • (2-12) 

It acts on the direct product of the auxiliary space and the quantum space in the nth site, C 2 (E)V n , 
and satisfies the Yang-Baxter commutation relations involving a rational R— matrix PI El Ej- 
Taking the product of iV Lax operators in the auxiliary space, one obtains the monodromy 
matrix 

t w („) = m»)...M«) = (^] Ztt)' (2 ' 13) 

with the operators A N (u), D N {u) acting on the quantum space of the model Vn- They satisfy 
the following Yang-Baxter relations PP 

B N (u)B N (v) = B N (v)B N (u), 
{v-u + i)A N {v)B N {u) = (v -u)B N {u)A N {v) + i A N (u) B N (v) , (2.14) 
(v — u — i) D N (v) B N (u) = (v — u) B N (u) D N (v) — iD N (u)B N (v). 
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The quantum determinant of the monodromy matrix (j2.13j) is given by ^Hj 

det q T N (u) = A N {u) D N {u + i) - C N {u) B N {u + i) = (u + is) N (u + i- is) N . (2.15) 

It follows from the Yang-Baxter relations that the auxiliary transfer matrix t^iu), defined as a 
trace of the monodromy matrix over the auxiliary space 

t N {u) = tr T N (u) = A N (u) + D N (u) , (2.16) 

commutes with itself for different values of the spectral parameter, with the Hamiltonian of model 
f!2.3|) and the operator of the total spin 

\t N (u),t N (v)} = ft N {u),H N ] = \t N {u),S] =0. (2.17) 

Substituting ()2.13j) into ()2.16j) and taking into account the explicit form of the Lax operator 
(12.12)1 one finds that tjv(w) is a polynomial of degree N in the spectral parameter u with the 
operator valued coefficients 

t N (u) = 2u N + q 2 u N ~ 2 + ... + q N . (2.18) 

One deduces from ()2.17|) that the operators q 2 , q~N form a family of mutually commuting, 
SL(2) invariant integrals of motion. Together with the total momentum of the system p, their 
eigenvalues q = (q 2 , ?at) form a complete set of the quantum numbers specifying the solu- 
tions to the Schrodinger equation (|2.8|) . As a consequence, the spectral problem ()2.8|) can be 
reformulated as 

t N {u)^ q>p (z 1 ,...,z N ) = t N (u)'9! qa ,(z 1 ,...,Zif), (2.19) 

with tAr(-u) being an eigenvalue of the transfer matrix ()2.18|) . The Hamiltonian of the model, 
Eq. ()2.3|) . can be obtained in a similar manner from the fundamental transfer matrix, which is 
constructed analogously to (|2.16|) and (j2.13|) from the Lax operators acting on the direct product 
of two copies of the quantum space V (8) V |2| . 

As follows from their definition, Eqs. (|2.18|) . (|2.16|) and (|2.12J) . the integrals of motion q^ take 
a form of polynomials in the spin operators S n (with n = 1, N) of degree k. For instance, 

q 2 = -2 J2 0rn ■ S n ) = -S 2 + Ns(s-1), (2.20) 

m>n 

with S 2 being the Casimir operator corresponding to the total spin of iV particles. In the 
representation (|2.5|) . is given by a complicated kth order differential operator and Eq. (|2.19j) 
leads to the system of (N — 1)— differential equations (g^ — qk)^ q , p (zi, . . . , z N ) = with k = 
2, . . . , N. Its exact solution for arbitrary iV becomes problematic. To overcome this difficulty we 
shall apply the method of Separated Variables (SoV) developed by Sklyanin in [TT| . 



3. Separation of Variables for the 51/(2, R) magnet 

Let us construct an integral representation for the eigenfunctions of the quantum SL(2,M.) mag- 
net, (z±, ... , zn), by going over from the coordinate z— representation to the representation 
of the Separated Variables (p, x) = (p,Xi, ...,Xjv_i) 

^ qtP (z 1 ,...,z N )= d N ~ 1 xii(x)U Pj:B (zi,...,z N )^ q (x), (3.1) 



5 



where /x(a?) is the integration measure on the a?— space and integration region will be specified 
below. U PjX (zi, . . . , zn) is the kernel of the unitary operator corresponding to this transformation 

U p , x (zi, . . . , z N ) = (zi, . . . , z N \p, x) , (3.2) 

where we introduced standard notations for the bra- and ket-vectors on the quantum space of 
the model. & q (x) is the eigenfunction in the separated coordinates 

5{p - p') = (p', x\q q , p ) = J V N z (U p , >x ( Zl , . . . , z N ))*^ q>p {zi, ...,z N ). (3.3) 

A unique feature of the SoV representation is that the eigenfunction & q {x) is factorized into a 
product of functions depending on a single variable & q (x) ~ Q(xi) . . . Q(xn-i)- 

In this Section, we shall obtain an explicit expression for the transition function to the SoV 
transformation, U PjX (zi, . . . , z^). The integration measure fi(x) and the properties of the eigen- 
functions & q (x) will be discussed in Sections 4 and 5, respectively. 



3.1. Basics of the SoV method 

To construct the unitary transformation to the SoV representation, Eq. ()3.2|) . we apply the 
Sklyanin's approach [H]. In this approach, the basis vectors \p, x) entering ()3.2j) are defined 
as eigenstates of the operator B^{u). We recall that B^(u) was introduced in (j2.13|) as an 
off-diagonal component of the monodromy operator. From the Yang-Baxter relations for the 
monodromy operator follows that [B N (u), B N {y )] = [B N (u), S-] = 0. This allows one to de- 
fine the eigenf unctions of the operator B^{u) in such a way that they are u— independent and 
diagonalize the operator of the total momentum (iS- — p)\p, x) = 0. 

One finds from (j2.13|) and (|2.12j) that Bjy(u) is a polynomial in u of degree N — 1 with 
operator valued coefficients, B^iu) = iS- tt^" 1 + . . .. It eigenvalues can be specified by the total 
momentum p and by the set of zeros x = (xi, . . . , xjv_i) 

B N (u) Up^zt, . . . , z N ) = p (u - - £at_i) U p , x {z u ...,z n ). (3.4) 

Notice that this relation is equivalent to the following system of equations 

U p>x (z) = p U PiX (z) , B N (x k ) U P:X (z) = 0, (k = 1, ...N - 1) , (3.5) 

where z = (zi, . . . , z N ). The solutions to Eqs. ()3.4|) and ()3.5|) are defined up to an overall 
normalization factor. It is convenient to choose it in such a way that U PtX (zx, . . . , zn-i) will be 
symmetric under permutations of any pair of the separated coordinates 

I 'p.-.r, .-■.,■„ •••(•- I • • • • , ZN-i) = U pr .. Xm ... Xn ...(zi, . . . , 2jV-l) • (3.6) 

Using ()2.12|) . (|2.13j) and (j2.11|) . one can verify that B^(u) is a self-adjoint operator on the Hilbert 
space of the model for real u 

(S J v(w)) t = B N {u*) . (3.7) 

This implies that the parameters x parameterizing its eigenvalues, Eq. ()3.4|) . take continuous 
real values. The set of corresponding eigenstates \p, x) is complete on the quantum space of the 
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model Vn and their orthogonality condition looks as 

(p',x'\p,x) = J V N zU PiX (zi, z N )(U pl;X '(zx, z N ))* 

= 8{p-p'){8(x-x') + .-.}-^^, (3.8) 

where 5(x — x') = Ylk=i $( x k ~ x 'k) an d ellipses denote the sum of terms with all possible 
permutations inside the set x = (xi, . . . , a?jv_i). 

The property ([3.7)1 also holds for the operators An(u), Dn(u), Cn(u) and for the transfer 
matrix, (tjv(w))t = t^{u*). Applying the both sides of the Yang-Baxter relations (|2.14jl to 
Up taB (zi, zn) and taking v = x k , one finds 

A N {x k ) U p>x (z) = a k (x) U PtX+iek (z) , D N (x k ) U p>x (z) = d k (x) U p>x „ iek (zi) . (3.9) 

Here the notation was introduced for the unit vectors e*, in the (N — 1)— dimensional x— space, 
such that (e fc ) n = 5 nk and x + ie k = (xi, . . . ,x k + i, . . . , x^v-i)- Substituting u = x k into ()2.15|) 
and applying ([3.9)1 . one finds that a k (x) d k (x + ie k ) = {x k + is) N (x k + i — is) N . The coefficients 
a k (x) and d k (x) depend on the normalization of U p ^ x (z), or equivalently on the definition of the 
integration measure in (|3.8jl . It is convenient to normalize U PjX (z) in such a way that 

a k (x) = A + (x k ) = (x k + is) N , 4(«) = A_(x fe ) = (x fc - is) N . (3.10) 

In Eq. (I3.9|l we have tacitly assumed that the function U PjX (z) can be analytically continued from 
real x into a finite strip in the complex plane. We will verify this property a posteriori. 

Using the solutions to ()3.4j) and (|3.5jl . one can decompose an arbitrary state on Vat over the 
basis of the functions U P;X (zi, . . . , z^). For the eigenstates of the model, the decomposition takes 
the form (j3.1|) and ()3.3|) . To obtain the wave function in the separated coordinates, $> q (x), one has 
to examine the action of the transfer matrix on the function U PtX (zi, . . . , zjy). According to ([2.18)1 . 
£jv(w) is a polynomial of degree N in u with the coefficient in front of u N_1 equal to zero. As a 
result, it can be reconstructed from its values at iV — 1 distinct points tjsr(xk) = A]y(x k ) +Djsr(xk), 
with k = 1, . . . , N — 1, using the Lagrange interpolation formula. Applying ()3.9[) one gets 

N-1 N-1 

t N (u)U PjX (z) = 2(u + ^x k ) Y]_{u - Xj)U p>x {z) 

k=l j=l 

N-1 

+ EII^^ A+(x k )U p , x+iek (z) + A4x k )U p (z) . (3.11) 
k=i j^k Xk ~ x i 

The wave function ([3.1)) has to diagonalize the transfer matrix, Eq. ([2.19)1 . Taking into account 
([3.11)1 . we find that the wave function in the separated coordinates satisfies a multi-dimensional 
Baxter equation 

t N (x k ) $ q (x) = (x k + is) N <£> q (x + ie k ) + (x k - is) N $ q (x - ie k ) , (3.12) 

with k — 1, . . . , N — 1. Here we assumed that the integration contour in ()3.1j) can be shifted 
into the complex x— plane. In addition, we took into account that the measure n(x) satisfies the 
following finite-difference equation (TTj 

H(x + ie k ) A + (x k ) yjXk-Xj + i 



fj,(x) A-(xk + i) ±1 x k -x 
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Figure 1: Diagrammatical representation of the function A u (zi, . . . , z^\w2, ■ ■ ■ , wjv)- The arrow 
with the index a that connects the points w and z stands for (z — w)~ a . 

with A±(x k ) defined in (J3.1(Jj) . This relation follows from the condition for £tv(m) to be a self- 
adjoint operator in the SoV representation, (p', x'\tjsr(u)\p, x)* = (p, x\tN(u*)\p' , x'). 
The same Baxter equation f)3. 12j) holds for a complex conjugated function 

%(x)8{p-p') = (* q , p \p',x). (3.14) 

To see this one uses f)2.19j) together with hermiticity of the transfer matrix 

t N (x k )(y qiP \p',x) = (t N (x k )ty q Jp' \x) = (y qiP \(A N (x k ) + D N (x k )\p',x) . (3.15) 

Here in the last relation A^{x k ) and D^{x k ) act on \p',x) as shift operators, Eq. (J3.9j) . Using 
(I3.14|) . one finds that $*(#) satisfies (|3.12|) and, therefore, one should expect that &* q (x) ~ &q( x ) 
for real x. We shall verify this relation below (see Eqs. ()5.17|) and (|5.8p ). 

As was mentioned in the Introduction, Eqs. ()3.13|) and ()3.12|) do not allow us to uniquely 
determine the integration measure and the eigenf unctions in the SoV representation since their 
solutions are defined up to a multiplication by an arbitrary periodic function f(x) = f(x±i). 
To fix the ambiguity one has to impose the additional conditions on the solutions to (|3.13|) and 
()3.12|) . As we will show in the next Section, these conditions follow from the explicit expression 
for the kernel of the unitary transformation to the SoV representation, U P)X (z). 

3.2. Construction of the SoV representation 

To determine the kernel U p ^ x (z) one has to solve the system (J3.5|) under the additional condition 
(|3.6|) . Constructing solution for U P;X (z), we shall follow the approach developed in Ref. [TJ] in 
application to the quantum SL(2, C) spin magnet. 

To begin with, one considers the first equation in the system (J3.5I) 

B N (x 1 )U P7X (zi, ■ ■ ■ , z N ) = . (3.16) 

Its solution can be found by making use of the invariance of the transfer matrix (J2.16)) under local 
gauge transformations of the Lax operators ()2.12|) and the monodromy matrix (|2.13|) [2UJ ITB*] 

L k (u) -»• L k (u) = Mfc- 1 L k (u) M k+1 , T N (u)^f N (u) = M{ 1 T N (u)M 1 . (3.17) 

Here M k are arbitrary 2x2 matrices, such that M^+i = M\ and detM^ ^ 0. Let us choose the 
matrices M k as 

Mk = ( J T 1 ) > M ^ = ( I 'f 1 ) < 3 - 18 ) 
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with wi, . . . , wn being arbitrary gauge parameters. It is straightforward to verify that being 
applied to the function 

(f) u (z k ; w k , w k+1 ) = (z k - w k )' s ~ iu (z k - w k+i y s+iu , (3.19) 

the gauge rotated Lax operator L k (u) takes the form of a lower triangle matrix 

L fc («) • ^(z fc ; Wfc+1 ) = ^ + {u _ %s)(j)u _ t j . (3.20) 

This suggests to define the following function 

N N 

Y u (z, w) = Yl <f>u{zk] Wk, Wk+i) = Y[( z k - w k )- 3 ' iu (z k - w k+i y s+iu . (3.21) 

k=l k=l 

Aside from the coordinates z — (z\, . . . , zn), it depends on auxiliary variables w = (wi, . . . ,Wn) 
and the spectral parameter u. Denoting the matrix elements of T N {u) as A N {u), . . . ,D N {u) in 
the same manner as ()2.13|) . one finds from ()3.20|) that Bn(u)Y u (z, w) = and 

A N (u)Y u (z, w) = (u + is) N Y u+i (z, w) , D N (u)Y u (z, w) = (u - is) N Y u _i(z, w) . (3.22) 

As follows from (|3.17|) . the operators An(u), . . . , Dn{u) are given by linear combinations of the 
operators A N (u), . . . , D N (u) with the coefficients depending on 1/wx. For w\ — » oo the two sets 
of the operators coincide since Mi = t and T^{u) = Tn(u). In this limit, the r.h.s. of (|3.21|) 
scales as Y u (z,w) ~ w\ s A u (z,w) with 

AJzi, . . 




H ( z » - w k )- s - tu (z k - w k+1 )~ s+m (z N - w N )~ s ~ iu . (3.23) 



k=2 



One deduces from 1)3.22)1 that for arbitrary W2, ■ ■ ■ , wn the function A u (z, w) satisfies the relations 

B N (u)A u (z, w) = 0, 

A N (u)A u (z, w) = (u + is) N A u+i (z, w) , (3.24) 
D N (u)A u (z, w) = (u- is) N A u _i(z, w) . 



It proves extremely useful to translate the obtained expression for the A— function, Eq. ()3.23j) . 
into a language of Feynman diagrams [Hj. Namely, one represents each factor [z — w)~ a in the 
r.h.s. of ()3.23j) by an arrow with the index a that starts in the point w and ends in z. In this 
way, the A— function can be represented as a Feynman diagram shown in Figure ^ 

Making use of the first relation in ()3.24|) , one finds that the general solution to Eq. ()3.16|) can 
be written as a convolution of the A— function with an arbitrary weight function Zn-i{w2, . . . ,wn), 
which is holomorphic in the upper half-plane and depends, in general, on the separated variables 



U PtX (z 1: . . . , z N ) = j V N 1 w A Xl (zi, . . . ,z N \w 2 , . . . ,w N ) Z N _ 1 (w 2 , . . . ,w N ) 

= [Ajv(zi) ® , (3.25) 



9 




Figure 2: Diagrammatical proof of the identity ()3.26|) . One inserts two vertical lines with the 
indices ±a = ±i(x 2 — X\) into one of the rhombuses and displaces them in the directions indicated 
by arrows with a help of identities shown in Figures [7| and |H1 

where D^ -1 ™ = Dto 2 . . . T>w^ with w n = u>*. Here the notation was introduced for the integral 
operator Aat(^i) with the kernel defined by the function A Xl (zi, . . . , zn\w 2 , ■ ■ ■ , wn)- 

Let us now require that the function U PiX (z) defined in ()3.25j) has to satisfy the second 
equation from the system ()3.5J1 . B^{x 2 ) U PjX (z) = 0. This leads to a rather complicated equation 
for the weight function Zn-i{w). To solve this equation, we shall propose a particular ansatz 
for Zjv-i(w) and verify that it leads to the expression for U PjX (z) which obeys (|3.5jl . Suppose 
that there exists such Zn-i(w) that the resulting expression for (|3.25|) is a symmetric function 
of Xi and x 2 . Then, the above requirement will be automatically fulfilled in virtue of ([3. 16)1 . The 
definition of Z^_\{yS) is based on the following remarkable identity 

[A N (xi) <g> A N - 1 (x 2 )](zi, . . .,z N \v 3 , ...,v N ) = [A N (x 2 ) <g> A N -x(xi)](zi, . . .,z N \v 3 , ...,v N ), 

(3.26) 

where the operator An(x 2 ) was defined for arbitrary A" in (13. 25 J) and (|3.23jl . The both sides of 
this relation can be rewritten as a convolution of two A— functions, Eq. 1)3.23)1 . In the Feynman 
diagram representation, the l.h.s. of ([3.26)1 is described by the diagram shown in Figure|21 Having 
the identity ()3.26|) in mind, we choose the function Zn-\(w) in the following form 

Z N -i(w 2 , ■ ■ -,w N ) = [Aat_i(x 2 ) <g> Z N ^ 2 ] (w 2 , . . .,w N ) , (3.27) 

where the function Z N _ 2 = Z N _ 2 (v 3 , . . . ,vn) is holomorphic in the upper half-plane, Imt> fc > 0, 
and depends, in general, on x 3 , . . . , x^. Substituting (|3.27|) into ()3.25j) and making use of ()3.26)) 
we find that U PtX (z) is invariant under permutations of x\ and x 2 

The proof of ()3.26j) can be performed diagrammatically without doing any calculation. It is 
based on elementary permutation identities shown diagrammatically in Figures [7| and |H1 Their 
derivation can be found in Appendix EI Similar identities have been also found for the Toda 
model [20] and the SL(2, C) spin chain ^3]. To prove (j3.26j) one inserts two auxiliary lines with 
the indices ±z(xi — x 2 ) into one of the rhombuses in Figure 121 Since the lines are attached to 
the same points and the sum of their indices equals zero, this transformation does not affect the 
l.h.s. of ()3.26j) . Then, one applies the permutation identity shown in Figure [7| and moves the line 
with the index i(x 2 — x{) to the left part of the diagram until it reaches the leftmost rhombus 
in which case one applies the identity shown in Figure |H1 Performing similar transformations on 
the second line with the indices — i{x 2 — x{), one moves it to the right part of the diagram. In 
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W N 




Figure 3: Diagrammatical representation of the function U x (z; wn)- The indices = s — ixk and 
Pk — s + iXk parameterize the corresponding factors entering 1)3.23)1 . The SL(2,M) integration 
()2.7j) over the position of internal vertices is implied. 

this way, one obtains the initial Feynman diagram but with the variables Xi and X2 interchanged, 
thus proving (|3.26|) . 

It is now straightforward to write a general expression for the function U PjX (z) satisfying the 
relations ()3.5j) and ()3.(i|) 



where z — (z\, . . . , zn) and U x (z; wn) is factorized into the product of N — 1 operators 



Here the operator k n (xN+i- n ) has an integral kernel A XN+1 _ n (zi, z n \w2, w n ) defined in 
f!3.23|) . It depends on a single separated coordinate x^ + i^ n (with n = 2,...,N) and being 
applied to an arbitrary function of W2, ■ ■ ■ , w n it increases the number of its argument by 1 (see 
Eq. ([3.25)) ) . The expression for the function U x (z; u>n), Eq. ()3.29|) . admits a simple diagrammatic 
representation. Replacing each A— operator in ()3.29j) by the corresponding Feynman diagram (see 
Figure HJ) one obtains that U x (z;wn) can be represented as the "Pyramide du Louvre" diagram 
shown in Figure El It consists of (N — 1)— rows with each row denoting a single A— operator in 
f|3.29|) . The k— th row of the pyramid consists of (N — A;)— lines carrying the index = s — ix^ 
and the same number of lines with the index /3). = s + ixk depending on the separated variable x^. 
Remarkably enough, a pyramid diagram similar to the one shown in Figure El describes the SoV 
transformation for the quantum SX(2, C) spin chain and the periodic Toda chain |13j . 

We would like to stress that a pyramid-like form of U x (z; wn) is a consequence of the factorized 
form of the kernel ()3.29|) . which is expected to be a general feature of the SL{2) magnets and 
related quantum integrable models. 




(3.28) 



U x (z; w N ) = [Ajv(zi) <8> A^-i^) <S> • • • <S> A 2 (x N - 1 )] (zi, z N \w N ) . 



(3.29) 
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Let us verify that the obtained expression for the function U x (z; wn), Eq. 1)3.29)) . verifies the 
defining relations ()3.5)) and l)3.fijl . Indeed, it follows from Eqs. (13.26)1 and (J3.16)) that U x (z;wn) 
is a completely symmetric function of x\, . . . ,Xj^-i, which is nullified by the operators B^{xk) 
for k — 1, . . . , N — 1. The integration over wn in (13.28)1 ensures that U PjX (z) is an eigenfunction 
of the operator of the total momentum iSL. To see this, one has to take into account that the 
kernel of the A— operator, Eq. ()3.23)1 . depends on differences of the coordinates and, therefore, 
the function U x (z] wn) is translation invariant, U x (z + e; wn + e) = U x (z; wn) with e real. 3 The 
normalization factor p Ns ~ 1 / 2 was inserted in the r.h.s. of ()3.28jl for the later convenience. 

By the construction, U x (z; wn) is a holomorphic function of z and w^- = w N in the upper-half 
plane. One finds from ()3.28)1 that U p ^ x (z) vanishes for p < since the integration contour over 
Re wn can be closed into the lower half-plane. Then, it follows from (13.1)1 that the eigenfunctions 
^ q, P (z) are different from zero only for p > 0. This is a common feature of states belonging 
to the Hilbert space ()2.9jl and (J2.tij) . The eigenstates with negative momenta, p < 0, can be 
constructed on the space of functions holomorphic in the lower-half plane. 

Since the kernel of the A— operator satisfies ()3.24jl . one verifies that, in agreement with (J3.9j) . 
the same relations hold for the function U p ^ x (z). Using ()2.16)1 one finds 

t N {xk) U P:X (z) = (x k + is) N U P:X+iek (z) + (x k - is) N U P;X _ iek (z) . (3.30) 

As was shown in Section 3.1, this leads to the mult i- dimensional Baxter equation for the wave 
function in the SoV representation, Eq. 1)3.12)1 . It follows from ()3.28)1 and ()3.29)1 that U PiX (z) is 
an entire function of x±, . . . ,xn-i and, as a consequence, the solutions to ()3.12jl have to possess 
the same property. Indeed, the analytical properties of U PiX (z) are in one-to-one correspondence 
with possible divergences of the Feynman integral in the r.h.s. of (J3.29)) . Using ()3.23jl one verifies 
that the integral is convergent for arbitrary complex x. 

3.3. Recurrence relations /Contour integral representation 

Expressions ()3.28)1 and ()3.29j) have a simple recursive form as a function of the number of particles 
N. Increasing this number, N — *■ N + 1, one has just to add an additional row to the pyramid 
diagram in Figure EH Namely 

U Xu ..., XN _ 1 (z 1 ,. . . ,z N ;w N ) = 

J V N ~ 1 v A XN ^ z N \v u . . . , U Xu ... :XN _ 2 (v 1 , . . . , wjv_i; w N ) , (3.31) 

with the A— function defined in ()3.23)) . The SL(2,M) integrals over the v— coordinates can be 
simplified by making use of the identity ()A.8)1 . In this way, one gets 

U Xl ,..., XN _ 1 {z 1 , ...,z N ; w N ) = [i 2s B(s + ix N -u s ~ ixN-i)] N+l 
N-i a 

x ] | / dr, rr^" 1 !! - rfe) 8 ^- 1 -^^..,^^^!), . . . , z^^r^!); % ) (3.32) 
k=i Jo 

where B(x, y) is the Euler beta-function and Zk{r k ) = (1 — r k )z k + r k z k+ i, so that z k (0) = z k and 
z k (l) = z k+ \. The [/—function in the r.h.s. is described by the same pyramid diagram but with 

3 Notice that the SX(2,R) integration in Eq. (|3.28|l can be replaced by the integration over the real axis 
J_ dwN. According to l|A.2l) . the result will differ by the normalization factor p 1 ~ 2s T(2s). 
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one row less. Its end-points are located in between the end-points of the pyramid in the l.h.s. of 
At iV = 2 one gets from (j3~23jl 

U Xl (z u z 2 - w 2 ) = (zi - w 2 )- s+iXl (z2 - w 2 )~ s - lXl ■ (3.33) 

Repeatedly applying ()3.32j) and using ()3.33j) as a boundary condition, one can express U x (z; wn) 
as a product of one-dimensional nested contour integrals. To save space, we do not present here 
the explicit expression. 

4. Integration measure 

To calculate the integration measure in the SoV representation, fi(x), one has to substitute the 
obtained expressions for the functions U PiX (zi, . . . , z^), Eq. ()3.28|) . into the orthogonality condi- 
tion ()3.8|) and perform integration. In spite of the fact that a mult i- dimensional integral in ()3.8|) 
seems to be rather complicated, it can be easily evaluated by making use of the diagrammatical 
representation for U PyX (z) (see Figure OJ). The analysis goes along the same lines as calculation 
of the integration measure for the SL(2, C) magnet [14]. 

Let us substitute (|3.28|) into ()3.8|) and consider the following scalar product 

(w' N , x'\w N , x) = J V N zU x (z 1 ,...,z N ]w N )(U x >(z 1 ,...,z N -,w' N ))*. (4.1) 

To evaluate ()3.8j) one has to perform a Fourier transformation of (w' N , x'\wn, x ) with respect to 
wn and w' N and multiply it by the additional factor p 2Ns ~ l , The function U x (z;wn) entering 
(14. 1|) is represented by the pyramid diagram shown in Figured To obtain (U x (z; w' N ))* one has to 
replace in Eqs. ()3.28|) and ()3.23|) the holomorphic "propagators" (z — w)~ a by complex conjugated 
expressions (z — w)~ a * = e lna * (w — z)~ a * . The function (U x (z; w' N ))* can be represented by the 
same pyramid diagram if one replaces in Figure El the indices otk and (3k by their conjugated 
expressions, ot* k = (3k and (3% — respectively, and flips the direction of all arrows. Each arrow 
is accompanied by the additional factor e J7rQ * or e J?r/3 fc. Combining these factors together along 
the rows of the pyramid diagram and using the identity + (3k = 2s, one finds that their total 
product equals e lnN ( N - 1 ) s _ As we will show below, in the final expression for the measure this 
factor cancels against a similar factor coming from the z— integration in (|4.1|) . 



w 2 




w' 2 



Figure 4: The scalar product of two pyramid diagrams at iV = 2. 
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It is convenient to flip horizontally the conjugated pyramid diagram, so that the point w' N 
will be located at the bottom of the diagram and the points z\, . . . , z^ at the top. The scalar 
product (|4.1|) is obtained by sewing together the pyramid and its conjugated counterpart at the 
points Zi, . . . ,Zn- The resulting Feynman diagram takes the form of a big rhombus built out of 
(N — l) 2 — elementary rhombuses (see the leftmost diagram in Figure 03 below) . Its tips have the 
coordinates wn and w' N . Notice that the indices a k and j3 k in the upper and lower part of this 
rhombus depend on two different sets of the separated coordinates x and x', respectively. 




(b) (c) (d) (e) 



Figure 5: Diagrammatical calculation of the scalar product of two pyramids. The fat points 
indicate the vertices which can be integrated out using the chain relation shown in Figure |H1 

Let us first calculate ()4.1|) at iV = 2. The corresponding rhombus diagram is shown in Figure^] 
Integration over Z\ and z 2 can be easily performed using the chain relation shown in Figure H3 
(see Appendix A for details) leading to (w 2 , x[\w2, ~ (w' 2 — ^2)°- Its Fourier transformation 
with respect to w 2 and w' 2 leads to the expression for (p',x[\p,Xi), which is divergent at p — p'. 
This means that, in agreement with (|3.8|) . (p',x[\p,Xi) should be understood as a distribution. 
To identify its form one has to regularize the corresponding Feynman integrals. As was shown 
in [Tl|, this can be achieved by shifting the indices as 



a k ->a k + e, f3 k -»• (3 k + e (4.2) 
and carefully examining the limit e — » 

(p',x'\p,x) = (pp') Ns - 1/2 lim J V N w N e ipWN J V N w' N e- ip '*» {w' N , x'\w N , x) e . (4.3) 

Repeating the calculation of the rhombus diagram in Figure EJ one gets 

(w' 2 ,x' 1 \w 2 ,x 1 ) e = e 2 ^^,^ ■ K - ^)- 4e , (4.4) 
where the notation was introduced for (see Eq. (|A.6j) ) 



(4.5) 



7 e (xi, x x ) = a(s — ix 1 + e, s + ix 1 + e) a(s + %X\ + e, s — ix 1 + e) 
q -2vks r ^ 2e + i( Xl _ x' 1 ))r(2e + i(x[ - xi))r 2 (2s) 
T(e + s + ixi)T(e + s — ixi)T(e + s + ix' 1 )T(e + s — ix[ 
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Then, one substitutes ()4.4j) into (|4.3)l . performs its Fourier transformation with a help of ()A.4[) 
and applies the identity 

km -i ^ 1 = 2n S (x) . (4.6) 

to get the expression for (p', x[ \p, x\) ~ S(p — p')5{x\ — x[), which matches the r.h.s. of Eq. ()3.8|) 
at N = 2 for 



ZlT 



r(s + ixx)T(s — ix x 



(4.7) 



r 2 (2s) 

This expression defines the integration measure in the SoV representation at iV = 2. It is 
interesting to note that /i(xx) coincides with the weight function for the continuous Hanh poly- 
nomials [211 . 

The calculation of the scalar product (|4.1|) for N > 3 is shown schematically in Figure El At 
the first step, one uses the chain relation, Figure^ to integrate out the left- and rightmost vertices 
with the coordinates z\ and z^, respectively. This transformation replaces two pairs of lines by 
two vertical lines with the indices ±i(xx — x[) (see Figure^) and brings the factor 7 e=0 (xi,x^) 
defined in (|4.5|) . Repeatedly applying the permutation identity, Figure we move one of the 
vertical lines horizontally through the diagram in the direction of another one until they meet and 
annihilate each other. In this way, we arrive at the diagram shown in Figure Et- In comparison 
with Figure E^i, it has two vertices less and the parameters x\ and x[ are interchanged. In this 
diagram there are already four vertices which can be integrated out using the chain relation, 
Figure |H1 Further simplification amounts to repetition of the steps just described. At each next 
step the number of vertices in the diagram is reduced and the a?— parameters got interchanged, 
x k *-* %j- Continuing this procedure, one obtains the diagram shown in Figure Eli, in which 
(AT— 1)— elementary rhombuses are aligned along the vertical axis. It is accompanied by the factor 
Y[j+k<N-i r Y e = ( x j> x 'k) defined in (j4.5p . Enumerating the rhombuses in Figure Eli starting from 
the top, one finds that the kth rhombus coincides with the N = 2 diagram shown in FigureE]upon 
replacement xx and x[ by x^-k and x' k , respectively. Based on the N = 2 calculations, one should 
expect that the chain of rhombuses produces the contribution ~ rifc=i $( x N-k — x' k ). This turns 
out to be correct, but in order to see it one has to regularize the Feynman integrals according to 
(14. 2|) and carefully examine their limit as e — > 0. For e ^ each elementary rhombus in Figure Eli 
can be replaced by a single line with the index 4e (see Figure Et)- This brings the additional 
factor n^i 1 le(%k, x N-k)- Finally, one integrates out the remaining (AT — 2)— vertices using the 
chain relation, Figure 01 and combines together all factors to find the following expression for 
the regularized scalar product (|4.1|) 

N-1 

(w' N ,x'\w N ,x) e = e ™ N{N - 1)s Y[ 7^0(^,4) Y[%(x k ,x' N _ k ) 

j+k<N-l k=l 

where A = 4e(N - 1) - 2s(N - 2). 

One substitutes (J4.8j) into (|4.3j) . performs its Fourier transformation with a help of (|A.4[) and 
applies (|4.6j) to get after some algebra 

JV-l 

(p',x>\p,x) = (2rr) N - 1 5(p-p')l[5(x k -x' N _ k ) 

k=i 
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wv2, .n. j<k r(i(*k - Xj))r(-i(x k - Xj)) 

x T JV (2s) Ar _ 1 — . (4.9) 

rijfc=i [r(s + ix fc )r(s - ixfc)] 

This relation agrees with Eq. (J3.8j) . but in distinction with the latter its r.h.s. is not symmetric 
in x. Indeed, simplifying the Feynman diagrams in Figures we have tacitly assumed that the 
7 e=0 — factors entering the r.h.s. of ()4.8|) are finite functions of x and x' . Taking into account ()4.5|) . 
one finds that this is true provided that Xj ^ x' k with j, k = 1, . . . , N — 1 and j + k < N — l. 4 
Under these conditions, only one term survives in the r.h.s. of ()3.8j) and it has the same form as 

Matching ()4.9j) into ()3.8j) one obtains the integration measure for N > 3 

N-1 N-1 

jj,{x) = c N Yl ( x k ~ xj) sinh(7r(x fc - Xj)) [ r (s + ix k )T(s - ix k )] N , (4.10) 

i,*=l k=l 

with the normalization factor c N = [T n2 (2s)(2tt) n - 1 (N - i)! 7r (A f -i)(^-2)/2j-i > The following 
comments are in order. 

In the SoV representation, the measure ()4.1U|) is a semi-positive definite function on the space 
of real as— variables. It vanishes at the hyperplanes Xj = x k . After analytical continuation to 
complex x, n{x) becomes a meromorphic function of x k (k — 1, . . . , N — 1) with the iVth order 
poles located along the imaginary axis at x k = ±i(s + n) with n 6 N. The measure decreases 
exponentially fast when x k goes to infinity along the real axis and the remaining x— variables 
take finite values 

V(x)~ e~ 2 ^\ x 2 k iNs ~ 1] , (4.11) 

as x k — > oo and lmx k = fixed. One verifies that the obtained expression for the measure, 
Eq. ()4.10|) . satisfies the functional relation ()3.13|) . 

The transition function U PjX (z), Eqs. (j3.28|) and (j3.29J) . satisfies the completeness condition 

dp / d N ^x fj,(x) (U PtX (wi, w N ))* Up^Z!, ...,z N )= K(z; w) , (4.12) 

where J d N ~ 1 x = Ylk=i dxj* and K.(z; w) is the reproducing kernel (kernel of the unity 
operator) on the quantum space of the model Vat (see Eq. (IA.7|) ) 

/ %TX S 
vNw n -( l - • • • > ^ = ■■■>**)> ( 4 - 13 ) 
k=1 [Zk w k) 

with ^(zx, . . . , zjy) holomorphic in the upper half-plane. The relation ()4.12j) is verified at N = 2 
by an explicit calculation in Appendix El For N > 3 the calculation is more involved and will 
be presented elsewhere. 



4 To restore the missing terms in the r.h.s. of Ij4.9|l . one has to use the symmetry of the pyramid in order to 
rearrange its rows and repeat the same calculation. 
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zi z-i zn zi z 2 z N -\ zn 

Figure 6: Reduction of the product of (N — 1) Baxter Q— operators to the pyramid diagram. 
The dotted lines shrink into a point at W\ — = . . . = wn~i- The dashed lines are reduced to 
a power of W\ in the limit W\ — > oo. 

5. Eigenfunctions in the SoV representation 

Let us consider the properties of the eigenfunction in the SoV representation <& q (x) defined in 
Eq. (jHHU). Using l|T3B|) . (jOJ) and ([X"2jl . one can rewrite (JHEJ) as 

(^,^1^) = y V N z (U x (z; w N )y^q,p(z) = e ipWN & q (x) ■ r(2g) ■ (5-1) 

As was shown in Section 3.1, the function $> q (x) satisfies the multi-dimensional Baxter equation 
()3.12|) . In this Section we shall construct the solutions to ()3.12|) . The analysis is based on 
the relation between the transition function to the SoV representation, U PiX (z), and the Baxter 
Q— operator for the SX(2,R) spin chain. 

By the definition [THj . the Q— operator acts on the quantum space of the model VV, depends 
on the spectral parameter u and satisfies the following conditions. It commutes with the trans- 
fer matrix ijv(^) and with itself for different values of the spectral parameter, [tjy(w), Q(f)] = 
[Q(w), Q(f )] = and fulfil the operator Baxter equation 

t N {u) Q(u) = {u + is) N Q{u + i) + (u- is) N Q{u - i) . (5.2) 

The eigenstates of the model ()2.8j) diagonalize the Q— operator 

Q(u)^f Piq (zi, ...,z N ) = Q q {u)^p^ q (zi, ...,z N ), (5.3) 

and the corresponding eigenvalues Q q {u) satisfy the same equation (|5.2|) . The Schrodinger equa- 
tion (|2.8|) is equivalent to (|5.3|) and the energy E q can be calculated as a logarithmic derivative 
of Q q {u) at u = ±is 122 D3|- 

The Baxter Q— operator for the homogenous SL(2,M) spin magnet was constructed as an 
integral operator in Ref. 



(u) (zi, ...,zn)=I V n wQ u (z 1 , ...,z n ;wi,.. .,w N ) • • • ,w N ) • (5.4) 
and the kernel was calculated following the Pasquier-Gaudin approach [20] as 

N 

Q u (z; w) = e™ sN J] (z k - w k+1 )- s+iu (z k - w k )~ s ~ lu , (5.5) 



fc=i 
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where wn+x = W\- The conjugated operator (Q(u*)y is defined in a similar way and its kernel 
is given by (Q u *(w ; z))* . Using (|5.5|) one can verify the following relation 

(Q(u*))t = PQ(u) = Q(n)P, (5.6) 

where P is the operator of cyclic permutations of iV particles, P^ = 1, 

P ^ q>p {zx, . . . , zjv-i, z N ) = ^q, p (z 2 , ...,z N , z x ) = e ldq ^q, p {z u . . . , z N _ lt z N ) . (5.7) 

Here in the second relation we took into account that the Hamiltonian of the model, Eq. ()2.3|) . 
is invariant under the cyclic permutations of particles and, therefore, its eigenstates possess a 
definite value of the quasimomentum 9 q = 2irk/N with k e N. It follows from ()5.6|) that the 
operators (Q(w*))t and Q(u) share the common set of the eigenf unctions and their eigenvalues 
satisfy the relation 

(Q q (u*))* = e i0 « Q q (u) . (5.8) 

At u = —is the kernel of the Q— operator ()5.5|) coincides with the kernel of the unity operator 
QUE} leading to 

Q(-is)=K, (5.9) 

so that Qq(-is) = 1. 

Notice that Eq. ()5.5j) looks similar to the definition of the kernel of the A— operator, Eq. ()3.23j) . 
which enters into the expression for the SoV transformation ()3.29j) . Namely, the expression for 
Q u {z ; w) e~ lnsN , Eq. (15. 5|) . coincides with the Y— function, Eq. (|3.21|) . and it differs from the 
A— function, Eq. ()3.23|) . by two factors, (z\ — Wi)~ s ~ lu and (zn — Wi)~ s+lu . This suggests that 
there should exist a relation between the transition function to the SoV representation U x (z; w^) 
and the kernel of the Q— operator. For the SL(2, C) magnet such relation has been found in 
Ref. dU. 

As a starting point, one considers the following transformation \I/ — > $ n , the so-called sepa- 
rating map [2"3] 

^ n ( Xl ,...,x N ^) = (a\Q(x 1 )...Q(x N _ 1 )\V) 

= J V N z\ ...J V N z N (n(^)YQ Xl (^ z 2 )... Q^izN-x; z N )^{z N ) , (5.10) 

where fl(zx, • • • , z^) is an arbitrary function holomorphic on the upper-half plane. When applied 
to the eigenfunction of the model, ^ Vl q(z\, . . . , Zji), it separates the variables for arbitrary Q(z) 

$ n (Xi, . . .,X N -i) = Q q (xx) . . . Qq(x N -x)(n\^q, P ) ■ (5.11) 

To reconstruct the eigenfunction ty qtP (zi, . . . , z^) one has to invert (|5. 1(J|) and define the inverse 
separating map <3? n — > Notice that for arbitrary Q(z) the transformation ()5.10p is not unitary 
and, therefore, it does not correspond to the SoV transformation ()5.1j) . Let us demonstrate 
however that there exists special state £l(z), for which the separating map ()5.10p becomes unitary. 
In that case, (J5.1(Jj) coincides with the SoV transformation (J5.1|) and the inverse separating map 
is defined by (|3.1j) . 

Following our diagrammatical approach, we use (J5.5)) and represent the kernel of the product of 
(N — 1) Baxter operators, [Q(xi) . . . Q(xjv-i)] (z; w) as the Feynman diagram shown in Figure 
to the left. Notice that it contains as a subgraph the pyramid diagram corresponding to the 
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transition function U PiX (z; w n) ■ It is possible to reduce the whole diagram to its subgraph as 
follows. One explores a freedom in choosing the w— coordinates to put w\ = W2 in the upper row 
of the left diagram in Figure |H1 The two lines connecting the points Wi and W2 with the vertex 
Hi merge into a single line with the index 2s, which corresponds (up to a numerical factor) to 
the reproducing kernel K(yi, Wi) = Q ms {y\ — wi)~ 2s , defined in Eqs. (jA.7|) and (j4.13J) . Therefore, 
one can effectively remove the y%— integration and put y\ = w\ instead. At the diagrammatical 
level this is equivalent to shrinking into a point the two lines connecting w\ and W2 with y\. In 
a similar manner, choosing w\ = . . . — wn-i in the upper row of the left diagram in Figure El 
one creates an avalanches of simplifications throughout the diagram which removes 2(k — l) lines 
in the kth row shown there by the dotted lines. The resulting diagram still contains 2(N — 1) 
additional lines (the dashed lines in the Figure EJ), which connect the point Wi with the vertices 
along the boundary of the pyramid. For W\ — > 00 their product scales as Wi N . To get rid 
of the additional lines one multiplies the whole diagram by w 2s ^ N ^ and sends w% — > 00. 
In this way, we arrive at the following relation 

U x (z; w N ) = i ~( N '^ N + 2 > lim w 2 ^-^ [Q( Xl ) . . . Q(x N ^)} (z; w 1} . . . , w 1} w N ) . (5.12) 

wi—*oo 

It can be rewritten in an operator form by introducing the special state 

N—l ins „7,2s p ins 

n m ,* N ( Zl , ...,z N )= i-V-WW J] 6 W ° ■ 2s , (5.13) 

^ (z k - w y s {z N - w N y s 

which depends on two complex parameters Wq and w^. This state is normalizable with respect 
to the SL(2, R) scalar product 

(^,^|^o,^) = ^(wq- 1 ~ ^" 1 )" 2s(JV " 1) « " ^)" 2s (5-14) 

and, therefore, it belongs to the quantum space of the model. Then, one finds from ()5.12j) the 
following relation between the unitary transformation to the SoV representation and the product 
of (N — 1) Baxter Q— operators 

U x {z;w N )= lim (z\ Q(xi) . . .Q(x N -i)\^w ,iD N ) , (5.15) 

wo— >oo 

where it is implied that one has to evaluate the matrix element and take the limit Wq — » 00 
afterwards. Substituting this relation into (|5.1j) and taking into account (|5.6j) . one obtains 

(w N ,x\S> qiP ) = lim (Q iIj0iiBN \Q(x 1 )...Q(x N - 1 )¥ N - 1 \^ q>p ) 

wo— > 00 

= c q (p)-j* w »Q q (x 1 )...Q q (x N - 1 ), (5.16) 

where c q (p) = e~ l9( i\im W0 _ >oo (Q € , ^ N=0 \^/ q ^). The p— dependence of c q (p) follows from (|5.1|) as 
c q(p) = c q ■ p _s ^ Ar ~ 2 - )_1//2 /r(2s). Matching the second relation in ()5.16|) into the r.h.s. of (|5.1|) one 
obtains that the wave function in the Separated Variables is given by the product of (N — 1) 
eigenvalues of the Baxter operator (J5.3j) evaluated for real u = x k 

$ q (x) = C q Q q (Xl) . . . Qq(x N -l) , (5.17) 

with the factor c q depending on the normalization of the eigenstates ^ q , p - 
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It is well-known [I""", [""| E2j that the eigenvalues of the Baxter operator for the quantum 
SX(2, R) magnet, Q q (u), are polynomials in u of degree h G N defined by the total spin of the 
system 

S 2 * q>p (z) = (h + Ns)(h + Ns - 1) tf, iP (z) , (5.18) 

or equivalently — —(h + Ns)(h+Ns — l) — Ns(s — l) in Eq. (|2.20|) . Then, one finds from (|5.17|) 
that the wave function & q (x) is given by the product of (N — 1) polynomials in the separated 
variables. The eigenfunctions of the model $ q , p (z) are orthogonal to each other with respect to 
the SL(2,M.) scalar product (|2.9jl . Going over to the SoV representation ()3.1|) and using ()3.8j) . 
one can write the orthogonality condition as 

(*Vj/I* w > = (^'I^sov^P-p') = *(p-lO<W , (5.19) 

where we took into account that the spectrum of the integrals of motion for the SL(2, IR) magnet is 
discrete [""2*1 "o"|. Here the notation was introduced for the scalar product in the SoV representation 

($ q '|$q) SoV = / d?- x xn(x)($j(x 1 ,...,x N - 1 )*$ q (x 1 ,...,x N -i). (5.20) 

Together with ()5.17j) and ()5.8j) this leads to the orthogonality condition on the space of solutions 
to the Baxter equation ()5.2|) and (|5.3|) 

» N-l 

/ d N - l xfi(x) TT Q q (x k ) Q q ,{x k ) ~ . (5.21) 

At = 2 this condition alone allows us to obtain the solution to the Baxter equation. Given 
that the integration measure at iV = 2, Eq. (|4.7jl . coincides with the weight function for the 
continuous Hanh orthogonal polynomials [21], one finds that Q q (x) equals the same polynomial 

Q»W = * F *{ 2s ,2s V' (5 - 22) 

with nonnegative integer h defined in ()5.18|) . This result is in agreement with the previous 
calculations [3 1221 UH IS]- For higher iV the solutions to the Baxter equation ()5.2j) can be 
expanded over the N = 2 solutions ()5.22|) with the coefficients satisfying the iV— term recurrence 
relations [221 013 • 

Substituting (|5.17|) into ()3.1|) . one obtains the integral representation for the eigenfunctions of 
the model ty qtP (z). As shown in Appendix [BJ the resulting expression for ty qp {z) coincides with 
the well-known highest-weight representation for the eigenstates in the Algebraic Bethe Ansatz 
approach "TJ, (J2J 

V q>p {z) =B N {\ 1 )...B N {\ h )n p {z). (5.23) 
Here the operator B^{u) is the off-diagonal element of the monodromy matrix ()2.13|) and the 



state ^l p (z) is defined as 



2s-l e iirNs f N 

Q p (z) = P v{2s) I ^ Vw e*» H(z k - w)~ 2s . (5.24) 



lmui>0 
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The parameters Ai, . . . , A^ are simple roots of the eigenvalues of the Baxter operator, Q q (u) = 
Y[k=i( u ~ ^k)- They satisfy the Bethe equations 



which follow from the Baxter equation (j5.2j) for Q q (u). The fact that the two different representa- 
tions for the eigenstates, Eqs. (j3.1|) and ()5.23|) . coincide implies that the Algebraic Bethe Ansatz 
is complete. We recall that the quantum space of the SL(2, M) magnet is infinite-dimensional for 
an arbitrary number of sites N and, as a consequence, the well-known combinatorial complete- 
ness analysis of the Bethe states (see Ref. [21] and references therein) is not applicable in that 
case. 

6. Conclusions 

In this paper we have studied the spectral problem for the quantum SL(2, R) spin magnet within 
the framework of the Quantum Inverse Scattering Method. This model has previously emerged 
in high-energy QCD as describing the spectrum of anomalous dimensions of composite high-twist 
operators. 

The central point of our analysis was the construction of the representation of the Separated 
Variables for the SL(2,~R) spin magnet. Following the Sklyanin's approach, we presented a gen- 
eral method for constructing the unitary transformation to the SoV representation. It allowed us 
to obtain the integral representation for the eigenfunctions of the model and calculate explicitly 
the integration measure defining the scalar product in the SoV representation. We demonstrated 
that the language of Feynman diagrams becomes extremely useful in establishing various prop- 
erties of the model. In particular, we found that the kernel of the unitary transformation to the 
SoV representation takes the factorized form (|3.29|) . In terms of Feynman diagrams this implies 
that the kernel can be described by the pyramid diagram shown in Figure H3 Notice that the 
same diagram has been previously encountered in the construction of the SoV representation for 
the SL(2, C) spin magnet in Ref. [T3J[T3|. The approach described in this paper can be applied to 
other quantum integrable models like the periodic Toda chain [2U1II3] and the DST model |25j . 
which represent degenerated cases of the «SX(2,IR) spin chain. 

Remarkably enough, many nontrivial properties of the SoV representation can be deduced 
in our approach from a few elementary identities between Feynman diagrams like the chain 
and permutation identities shown in Figures El and [7J respectively. We found that the kernel 
of the unitary transformation to the SoV representation is given by the product of the Baxter 
Q— operators projected onto a special reference state. This allowed us to establish the relation 
between the wave functions in the separated variables and the eigenvalues of the Q— operator 
for the SL(2, M) spin magnet. Using the well-known fact that the latter are polynomials in the 
spectral parameter, we demonstrated the equivalence between two different expressions for the 
eigenstates obtained within the Algebraic Bethe Ansatz and the SoV method. 

So far we have assumed that the spin chain was homogeneous. Another advantage of our 
method is that it can be easily extended to the case of the inhomogeneous closed spin chain. 
For the SL(2, M) chain with the spin in the kth site equals Sk, the unitary transformation to 
the SoV representation is defined by the same pyramid diagram (see Figure EJ) with the a— and 




(5.25) 
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Figure 7: Permutation identity. 




Figure 8: Special case of the permutation identity. It is obtained from Figured by sending one 
of the external points to infinity. 

(3— indices modified in the following way. For Sk ^ Sj, the indices take different values for the 
different lines in the same row. If one denotes by a n k and (3 n k the indices carried by the fcth pair 
of adjacent lines in the nth row from the bottom, (1 < n < N — 1 and 1 < k < N — n), then 
&nk = Sk — ix n and (3 n k = s n+ k + ix n . In addition, in order to preserve the SL(2, R) invariance of 
integrals, one has to modify the integration measure ()2.7|) in all vertices of the diagram. Namely, 
one has to replace (2 \mw n k) 2s ~ 2 — > (2lmw n k) ank+l3nk ~ 2 in the integration measure (J2.7)) at the 
vertex w n k, to which the lines a n k and (3 n k are attached. 

Acknowledgements 

This work was supported in part by the grant 00-01-005-00 of the Russian Foundation for Fun- 
damental Research (A.M. and S.D.), by the Sofya Kovalevskaya programme of Alexander von 
Humboldt Foundation (A.M.) and by the NATO Fellowship (A.M.). 

A Appendix: Feynman diagram technique 

In this Appendix we collect some useful formulae, which were used to prove the permutation 
identity shown in Figured and to calculate of the integration measure (j4.10|) . 
To evaluate the 5L(2,R) integrals we use the identity 




(A.l) 
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with T>w defined in (|2.7j) . w = w* and p > 0. The momentum representation for function ^(z) 
holomorphic in the upper-half plane is defined as 

JO 1 Jlmz>0 

Applying the integral representation for the propagators 

i p—ina/2 roc 

dp e i P (z-w) p a-l ) ( A 3 ) 



r(a) 

which is valid for lm(z — w) > 0, one finds 



/ Vw - — = 9(p) p a - 2s • e— / 2 -LI , (A.4) 

Jim W >o {z-w) a T{a) 

Vw{z -w)~ a {w -v)- 13 = a(a,l3)(z-v)- a - p+2s . (A.5) 

lmu>>0 

Here the notation was introduced for 

_ r(a + /3-2 g )r( 2g ) 

e r(a)rcs) • ( 6) 

Eq. (|A.5p is the "chain relation" shown in Figure El From Eq. (jA.4j) and ()A.2|) one finds that 

v w -^—n W ) = nz). (A.?) 



Using the Feynman parameterization 

x 1 "x 2 " = ^ttvfT^ ^ drr°- 1 (l-r) /J - :L [rxi + (l-r)x 2 ]- a ^ 



^ _ r(a + /?) /- 1 
r(a)r(/3) 



and applying (|A.7|) . one finds that for a + (3 = 2s [TH] 

/ Vw- JlZ y( w) = r ( 2a ) / 1 dTr a-i (1 _ T) ^ (TZi + (1 _ r)z2) _ (A8) 

iim«,>o \zi ~ w) a (z 2 - w)P T(a)T((3) J 

The permutation relation shown in the Figure [7| involves the Feynman integral 

I(z, v;xi,x 2 ) — [ Vw- r — r^—. -5— r — , (A. 9) 

Ji mw >o (w-v^tw-v^^-wyi^-wY^ 

with a fc = s — ixk and fik = s + ix^. One applies (|A.8|) for ty(w) — (w — Vi)~ ai {w — v 2 y^ and 
calculates the integral in terms of 2 Fi — hypergeometric function as 

I(z, v ; Xl , x 2 ) = e' ms (z 2 - v 2 )~ a > ( Zl - v x y^ (z, - v 2 )^ +a *~ 2s 2 F X (a u a 2 ; 2s; £) , (A.10) 

with £ = [z\ — z 2 )(vi — v 2 )/[(zi — Vi)(z 2 — v 2 )}. From here, one gets for X\ — x and x 2 = y 

{z 2 - v 2 y^I(z, v;x,y) = ( Zl - J(z, v;y,x). (A.ll) 
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Figure 9: Chain relation. 



Multiplying the both sides of Eq. (jA.llj) by v 2 ~ lx and taking the limit v 2 — > oo, one obtains the 
reduced permutation relation shown in Figure |H| 

Finally, let us verify the completeness conditions ()4.12j) at iV = 2. We use ()A.3j) to rewrite 
the transition function U PiX (zi,Z2) defined in ()3.33j) and ()3.28j) in the momentum representation 

as 

^, M= i^r* r^(?i)\~ S(Vl+n - p) . (A , 2) 

Y{s-ix)Y{s + ix) J p\ J p 2 \PiJ 

At N = 2 the pyramid diagram in Figure E3 consists of two lines. The variables pi and p 2 have the 
meaning of the momenta that flow along these lines. Using similar expression for the conjugated 
function (U P)X (wi,w 2 ))* as a double integral over the momenta p[ and p' 2 , one examines the 
product n(x)U PtX (zi, z 2 ){U p ^ x (w\, w 2 ))* with the integration measure /i(x) given by ([4.7)1 . One 
notices that the x— dependence of the integration measure is compensated by the T— functions 
entering the r.h.s. of (|A.12|) . As a result, the p— and x— integration in (|4.12|) can be easily 
performed leading to the following combination of the 5— functions 

<Kpi+P2 -p[ -p' 2 )s(ln?^j ) = PlP2 S {p[ - pi)S(p' 2 - pa) ■ (A.13) 

V P2P1 J P1 + P2 

Combining together all factors one evaluates the l.h.s. of (j4.12j) as 

f°° (r> r> ^ 2s ~ 1 _ _ p i,rs p i,rs 

dpi / dp 2 [PlP 2 2) ^-^M^) = . (A.14) 

Jo T 2 (2s) (zj - w^ 23 (z 2 - w 2 ) 2s 

This expression coincides with the r.h.s. of ()3.28|) at = 2. 



B Appendix: Relation to the Algebraic Bethe Ansatz 

In this Appendix we demonstrate an equivalence between two different representations for the 
eigenstates of the SX(2, IR) magnet, Eqs. (|3.1j) and (|5.23|) . obtained within the SoV method and 
the Algebraic Bethe Ansatz, respectively. 

To begin with, one constructs the state equal to the product of the highest weights in sites 

N 

n(z) = l[e™zf s , S+n(z) = 0. (B.l) 

k=l 

It is analogous to a pseudovacuum state for compact spin magnets. This state annihilates the 
lower off-diagonal matrix element of the Lax operator (|2.12p and, therefore, diagonalizes the 
transfer matrix (j2.16j) 

t(u)Q(z) = [(u + is) N + (u- is) N ] Q(z) . (B.2) 

Due to the SL{2) invariance of the transfer matrix, \t(u), S] = with S = J2k=i $k being the total 
spin of the magnet, the state e wS Q(z) = fl(zi — w, . . . , z N — w) also verifies ()B.2|) 5 . Applying 
5 Notice that in distinction with (|B.ip . the state Sl(z\ — w, . . . , z» — w) has a finite SL(2,M.) norm for Imw < 0. 
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(|A.2J) . one performs its Fourier transformation 



n 2s-l 

n p {z) = ^— - I Vw e ipw Q( Zl - w, . . . , z N - w) . (B.3) 



r(2a) 



Im w>0 



The state Q. p {z) carries the total momentum p > 0, diagonalizes the transfer matrix (jB.2|) and 
satisfies the following relations 

is- n p (z) =pn p (z) , s 2 n p (z) = s( s - i)a p (z) . (ba) 

Its SX(2,R) norm can be easily calculated from ()A.4|) and (|A.1|) as 

(n p ,\a p ) = J v N z(n p ,(z)yn p (z) = 6(p-p')^——. (b.s) 

This implies that the state fl p (z) belongs to the quantum space of the model. Then, it follows from 
Eqs. (IB. 2 J) and (|B.4J) that, in virtue of complete integrability of the model, Q p (z) is an eigenstate 
of the Hamiltonian (|2.3|) . Its total SL(2, R) spin can be found from (jB.4|) and (|5.18|) to be h = 0. 
According to ()2.4j) and (|2.3p . the energy of this state is E q = since (J n ,n+i — 2s)fl p (z) = for 
n = 1, . . . , N. The values of the corresponding integrals of motion q can be found by matching 

([B~2|i into iprraji . 

Let us transform the state f2 p (z) into the SoV representation. One finds from (|5.1fij) 

(w N ,x\a p )= lim (Simoon \ Q{x x )... Q(arjv_i)|n p ), (B.6) 

since Pf2 p (i*) = £] p (z) (see ()5.7p and (|B.ljl ). Notice that fi(zi — w, . . . , z N — w) is proportional to 
the reproducing kernel ()4.13j) and its convolution with the kernel of the Q— operator, Eqs. ()5.4j) 
and (|5.5jl . can be easily performed leading to Q(u)\Q p ) = \Q P ). In this way, one finds from ()B.6j) 
and 

1 (2s) 

Together with the completeness condition (j4.12j) this leads to 

O (g) = ^.-1/2^^-1)^44). /" ^-Ig. ^ x )U ptW (g) . (B.S) 



Substituting this relation into (jB.5|) and taking into account ()3.8|) . one gets 

d^cc u(sc) = - 1 - , (B.9) 

where the measure n(x) is given by (|4.10|) . 

Let us apply the operator B N {Xi) . . .B N (X h ) to the both sides of ()B.8|) . with B N being the 
off-diagonal matrix element of the monodromy operator (j2.13|) and some parameters. By the 
definition (|3.4j) . the transition function U PtX (z) diagonalizes the operator 

N-l ( h ~\ 

fltf(Ai) . . . B N (\ h )U PjX (z) = (-l)Miv-iy "Q J _ Xk) I Upx ^ (B 1Q) 

j=l U=l J 
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Let us choose Ai, . . . , Xh to be roots of the eigenvalue of the Baxter Q— operator, Q q (\k) = 0, so 
that the expression in the parenthesis coincides with Q q (xj). Finally, one gets from ()B.8j) and 




(B.ll) 



with the normalization factor c(p) = p^«+/»-i/2(_i)A(JV-i)j(iV-i)(Jv+4) a- The l h s of flgjj]) coin _ 
cides with the well-known expression for the eigenstates in the Algebraic Bethe Ansatz, whereas 
the r.h.s. defines the integral representation for the same eigenstate in the SoV representation, 
Eq. (EU). 
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